Triangular flow in event-by-event ideal hydrodynamics in Au+Au collisions at 

v /5^ = 200A GeV 
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The first calculation of triangular flow 113 in Au+Au collisions at ^/snn = 200^4 GeV from an 
event-by-event (3+l)-d transport+hydrodynamics hybrid approach is presented. As a response to 
the initial triangularity €3 of the collision zone, ^3 is computed in a similar way to the standard 
event-plane analysis for elliptic flow V2. It is found that the triangular flow exhibits weak centrality 
dependence and is roughly equal to elliptic flow in most central collisions. We also explore the 
transverse momentum and rapidity dependence of v 2 and V3 for charged particles as well as identi- 
fied particles. We conclude that an event-by-event treatment of the ideal hydrodynamic evolution 
starting with realistic initial conditions generates the main features expected for triangular flow. 

PACS numbers: 25.75.-q,24.10.Lx,24.10.Nz 



Collective behaviour of particles emitted from relativis- 
tic heavy ion collisions, such as elliptic flow, is one of the 
earliest predicted observables that indicates fluid-like be- 
haviour of the created hot and dense nuclear matter [lj- 
4]. The pressure gradients need to be large enough to 
translate an early stage coordinate space asymmetry to 
a final state momentum space anisotropy. Therefore, the 
high values of the second coefficient of the Fourier ex- 
pansion of the azimuthal distribution of the particles, i>2, 
that have been reported by the experiments at the Rela- 
tivistic Heavy Ion Collider (RHIC) [Hi] have led to the 
conclusion that the quark gluon plasma is a nearly per- 
fect liquid [!G3. 

The collective flow observables manifest themselves in 
multi-particle-correlations as well. For example, the so 
called cumulant method has been very successful in quan- 
tifying the harmonic coefficients of the azimuthal parti- 
cle distributions [llj. Recently Arj-Acj) correlations have 
been explored in a new manner by extracting a triangular 
flow signal from the data that is responsible for most of 
the structures that were previously attributed to other 
mechanisms [T3, EH- Features like long-range rapidity 
correlations on the near- and away-side accompanied by 
a conical structure on the away-side have been often re- 
ferred to in the context of jet-medium interaction. The 
preliminary PHOBOS data show a long range correlation 
in rapidity which would be supported by an initial state 
generated from a flux tube picture like in NEXspheRIO 
jlil [l5j . In this model one is also able to observe the 
features like the ridge and the "cone" in the two-particle 
correlations. 

The triangular flow, V3, is the third coefficient of the 
Fourier expansion of the azimuthal distribution of the 
final particles in momentum space with respect to the 
corresponding event plane angle ^3 that is defined be- 
low. This angle fluctuates randomly event by event in 
contrast to the well-known event plane angle ^2 used for 
the elliptic flow analysis, which is strongly correlated to 
the reaction plane. The triangular flow is assumed to 
be the response to a triangular shape of the initial state, 
described by a nonzero triangularity, £3 , that arises from 



the fluctuations of the initial collisions. 

In contrast to the extensively studied observables like 
directed flow (vi), elliptic flow (i^) and the fourth har- 
monic (1*4), the triangular flow (1)3) is not correlated to 
the reaction plane that is defined by the beam axis and 
the impact parameter axis of the collision. The initial 
state fluctuations in the transverse plane are random with 
respect to the reaction plane. In a standard hydrody- 
namic calculation with smooth initial conditions only the 
even coefficients of the Fourier expansion are non-zero at 
midrapidity. The odd coefficients vanish by symmetry in 
collisions of identical nuclei which is the reason why they 
have not been studied so far. 

In this paper the latest version of the Ultra-relativistic 
Quantum Molecular Dynamics (UrQMD) [H, [l7| to- 
gether with ideal relativistic fluid dynamics is used to 
explore this new observable [32| . The full event- by- event 
setup of this hybrid approach allows to extract a V3 com- 
ponent from the final particle distributions. The method 
is very similar to the standard elliptic flow event plane 
measurement and will be outlined below. Predictions 
for the impact parameter dependence and the transverse 
momentum dependence of identified particles are made. 

Let us now review the main ingredients of the hybrid 
approach [l8|, that are relevant for the development of 
triangular flow. The initial binary nucleon-nucleon colli- 
sions are modeled in UrQMD following the Lund model of 
nucleon-nucleon reactions [20J involving color flux tubes 
excitation and fragmentation processes that provide long 
range rapidity correlations and fluctuations in the energy 
deposition in the transverse plane. For Au+Au collisions 
at the highest RHIC energies the starting time for the hy- 
drodynamic evolution has been chosen to be £ s tart = 0.5 
fm in order to fit the final pion multiplicity at midra- 
pidity. Only the matter around midrapidity (|y| < 2) 
is considered to be locally thermalized and takes part in 
the ideal hydrodynamic evolution. The more dilute spec- 
tator/corona regions are treated in the hadronic cascade 
approach throughout the reaction. To map the point par- 
ticles from the UrQMD initial state to energy, momen- 
tum and net baryon density distributions each particle is 
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represented by a three-dimensional Gaussian distribution 

Ei- 

The ideal hydrodynamic evolution [22|, for the hot 
and dense stage of the collision translates the initial fluc- 
tuations in the transverse energy density to momentum 
space distributions. A hadron gas equation of state [24| 
has been used because we are aiming here only at qual- 
itative statements and not at quantitative comparisons. 
This equation of state (EoS) has been extensively tested 
and gives reasonable results for multiplicities and parti- 
cle spectra. Furthermore, given the same average speed 
of sound during the evolution, flow observables are not 
sensitive to the details of the EoS [IE HE) ■ 

The transition from the hydrodynamic evolution to the 
transport approach when the matter is diluted in the late 
stage is treated as a gradual transition on an approx- 
imated constant proper time hyper-surface (see [H, [27j 
for details). The final rescatterings and resonance decays 
are taken into account in the hadronic cascade. 

The above event-by-event setup includes all the main 
ingredients that are supposed to be necessary for the 
build up of triangular flow. Since the complete final state 
particle distributions are calculated, an analysis similar 
to those applied by experimentalists is used. 

The definition of the participant eccentricity can be 
generalized to the triangularity defined as 



yj (r n cos(ruj))} 2 + (r n sm(n^))' 



(1) 



where in contrast to [12| the factor in front of cos(3(/>) is 
taken to be r 3 since this is the more consistent way to do 
the Fourier expansion of the distribution function, (r, </>) 
in this case are polar coordinates corresponding to the 
transverse plane in coordinate space. We have calculated 
the triangularity of the UrQMD initial state for all the 
particles that are newly produced or have undergone at 
least one interaction and with a rapidity between — 2 < 
y < 2 at the starting time of the hydrodynamic evolution 
(tstart = 0.5 fm). The centrality dependence and the 
values of the eccentricity as well as the triangularity are 
comparable to those obtained by a Glauber Monte Carlo 
approach [l3[ (see Fig. @|. 

The particle distribution in coordinate space is then 
transferred to the final state particle distribution in mo- 
mentum spaceby the pressure gradients during the hydro- 
dynamic evolution. Experiments are only able to mea- 
sure the momenta of the particles but not the coordi- 
nates, therefore, one has to find a way to generalize the 
elliptic flow analysis to triangular flow analysis in a con- 
sistent way. We propose here to use the standard event 
plane method [28[ and define an event plane for triangu- 
lar flow in the following way 



1 (pr sm(nd>„)) 
tf„ = -arctan-p K -^- 

n {pr cos(n(pp)) 



(2) 
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FIG. 1: (Color online) Distribution of event-plane angles ^2 
(dotted and full line) and $3 (open squares and circles) with 
respect to the reaction plane in central (6 < 3.4 fm) and mid- 
central (6 = 5 — 9 fm) Au+Au collisions at ^/snn = 200A 
GeV. 

As all the details of the final particles are known in our 
calculation, the resolution of the event plane angle can 
be improved by taking into account all the particles (also 
neutral particles) in a certain kinematic range (\r)\ < 2) to 
determine the event plane. The resolution is calculated 
by comparing two sub events and turns out to be around 
0.8 radians for the triangular flow plane angle and 0.7 to 
0.95 for the vi event plane angle depending on the cen- 
trality. This corresponds to an uncertainty in the angle 
of approximately 12 degrees in mid-central collisions. 

The distributions of the resulting event plane angles 
^2 and ^3 with respect to the known reaction plane 
(positive x-direction defines \V = 0) are shown in Fig. 
[T] As expected the elliptic flow event plane is Gaus- 
sian distributed and therefore correlated to the reaction 
plane, especially in less central events. The triangular 
flow plane shows a flat distribution between -60 and +60 
degrees since the fluctuations that lead to a u 3 component 
are random with respect to the reaction plane. These in- 
ternally consistent results increase our confidence in the 
analysis method proposed here. 

To explore the correlation between the initial spatial 
event plane and the corresponding final event plane we 
define an initial event plane angle $„ in analogy to Eq. 
©: 



1 (r™ sin(n0)) 
= — arctan- - — — - 

n (r n cos{n<p)) 



(3) 



where (j>t, i'p) are polar coordinates in momentum space. 



For elliptic flow this angle is the one that defines the 
so called participant plane. With the used conventions 
<&„ is defined in the region between —ir/n and +n/n. For 
convenience, we introduce a shifted angle $' n = + ir/n 
[mod 2ir/n], defined in the domain —ir/n < <f>^ < +tt/ti. 
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FIG. 2: (Color online) Correlation of the final event plane 
angles $2 (dotted and full line) and ^3 (open squares and 
circles) with the corresponding initial event plane angles $ 2 
and $3 in central (6 < 3.4 fm) and mid-central (6 = 5 — 9 fm) 
Au+Au collisions at ^/snn = 200 A GeV. 

The correlation between the final ($„) and initial ($' n ) 
event planes for two different centrality classes is shown 
in Fig. [2] There is a strong correlation in all four 
cases which has a similar shape. For elliptic flow there 
is a stronger contribution from the collision geometry 
that results in a centrality dependence of the correla- 
tion (^2 — $2) [29]- in more central collisions most of 
the elliptic flow comes from fluctuations while in more 
peripheral events the almond-shape geometry of the col- 
lision zone has a major effect on V2 which leads to an 
even stronger correlation of the two angles. For the tri- 
angular flow the distribution is very similar for central 
and mid-central events since it is only caused by fluctu- 
ations in both cases. This analysis has been carried out 
on an event by event basis and confirms that the final 
triangular flow is related to the initial triangularity. 

Further evidence for a strong correlation between 
initial state geometry and the final state momentum 
anisotropics is found by looking at the flow coefficients as 
a function of the corresponding e n . The flow coefficients 
are calculated by using the following formula 

v n = (cos{n((f> p - * n ))) (4) 

where it is important to note that the particle that is 
correlated to the event plane is removed from the event 
plane determination to eliminate auto-correlations. The 
final results for v n are obtained by dividing the above 
values by the event plane resolution of the corresponding 
centrality class. The same procedure has been applied in 
the event-by-event approach of Holopainen et al. 30] . 

In Fig. |2] we show ^2(62) and ^3(63) for two differ- 
ent centrality classes (b < 3.4 fm and 6 = 5 — 9 fm). 
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FIG. 3: (Color online) Dependence of v% and V3 on the initial 
£2 and £3 of charged particles in central (6 < 3.4 fm) and mid- 
central (6 = 5—9 fm) Au+Au collisions at ^/snn = 200^4 GeV 
at midrapidity (|^| < 1). 



All the curves behave linearly within error bars, a larger 
initial eccentricity /triangularity leads to a larger ellip- 
tic/triangular flow. For V3 the translation of the initial 
anisotropy into the final anisotropy is less efficient for 
non-central collisions than for central collisions. For el- 
liptic flow the lines have a similar slope, so the elliptic 
flow response is similar in central and more peripheral 
collisions. 

The impact parameter dependence of the flow coeffi- 
cients for charged particles is shown in Fig. |4] First of all, 
the triangular flow has a finite value which exhibits only 
a weak centrality dependence. This is another hint that 
V3 is only induced by fluctuations in contrast to i>2 which 
has a geometry influence in addition that is centrality 
dependent. 62 and €3 increase as a function of impact 
parameter since the almond shape of the overlap region 
gets more pronounced in peripheral collisions and the 
fluctuations are more important in smaller systems. For 
the third Fourier coefficient the ratio of ^3/63 decreases 
slightly for more peripheral collisions which reflects the 
shorter duration of the hydrodynamic evolution. There 
is less time to translate the initial state anisotropy to a 
final state momentum anisotropy. Furthermore, the rel- 
ative elliptic flow is larger than the relative triangular 
flow, so the transfer of coordinate space anisotropy to 
momentum space anisotropy is more efficient for lower 
harmonics than for higher ones. The same result has al- 
ready been found for Vi that is much smaller than V2 
31]. 

The momentum dependence of W3 and V2 for charged 
particles is shown in Fig. [SJ It shows that the elliptic 
flow is almost twice as large as the triangular flow in mid- 
central collisions. The elliptic flow results extracted with 
the event plane method are nonzero even in the most cen- 
tral collisions due to the fluctuations of the participant 
plane with respect to the reaction plane. In central colli- 



4 




0.1 
0.05 
0.0 



-* 

p a gj n a m am a a a a a m a a ° ° m " 



123456789 

b [fm] 

FIG. 4: (Color online) Impact parameter dependence of «3 
(open squares) and «2 (open circles) of charged particles in 
Au+Au collisions at ^/snn = 200A GeV at midrapidity (\r)\ < 
1). The full and the dotted line represent the ratios of V3/e3 
and V2 j to, respectively. 



sions the magnitude and the shape of the triangular flow 
is similar to the elliptic flow. This can be traced back 
to the fact the distribution of spatial event plane angles 
P(^2) and P(^3) are flat in central collisions (b=0 fm), 
rendering no correlation between vi jvz and the reaction 
plane. The remaining correlation for central collisions 
(see Fig. [1} of ^2 with the reaction plane arises because 
there is a contribution from finite impact parameter cal- 
culations (6 < 3.4 fm) in the most central class of events. 
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FIG. 5: (Color online) Transverse momentum dependence of 
V3 and V2 of charged particles in central (b < 3.4 fm) and mid- 
central (b = 5 — 9 fm) Au+Au collisions at -^/snn = 200^1 GeV 
at midrapidity (|^| < 1). 

The transverse momentum dependence of triangular 
flow for identified particles is shown in Fig. [6] We 
have checked that the elliptic flow for identified parti- 
cles is compatible with the experimental data. The fire- 



ball created at the highest RHIC energies is dominated 
by mesons and thus the pion flow is very similar to the 
charged particle flow. For the protons, the same mass 
splitting effect is seen for V3 as for elliptic flow (26[. In 
addition, the proton V3 is almost equal for central and 
mid-central collisions. 
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FIG. 6: (Color online) Transverse momentum dependence of 
V3 for pions and protons in central (b < 3.4 fm) and mid- 
central (6 = 5 — 9 fm) Au+Au collisions at -^snn = 200^4 
GeV at midrapidity (\rj\ < 1). 

Fig. [7] shows the pseudorapidity dependence of the 
two flow coefficients for charged particles in two central- 
ity classes. Due to the initial conditions generated by 
UrQMD with its reliance on fluxtube fragmentation one 
obtains a long-range A77 correlation that can be observed 
in the final state. The elliptic flow results are flat for at 
least two units of pseudorapidity whereas the triangular 
flow distribution is almost flat over the whole pseudora- 
pidity range (A77 = 4) covered by the present calculation. 

In conclusion, we have presented the first calculation 
of triangular flow from a (3+l)d ideal hydrodynamics 
approach in Au+Au collisions at */5nn = 200A GeV. 
The fluctuating initial conditions and the event- by- event 
setup are crucial for this observable. By translating ini- 
tial state triangularity to the final state momentum dis- 
tributions via pressure gradients a finite value of the third 
coefficient of the Fourier expansion of the azimuthal dis- 
tribution of the particles in the final state is generated. 

Our method is based on a generalization of the stan- 
dard event plane analysis that has been used for elliptic 
flow and can therefore be done in experiments in exactly 
the same way. While vi shows a strong impact param- 
eter dependence, i>3 exhibits only a weak centrality de- 
pendence and is close to i>2 in central collisions. The 
transverse momentum dependence of W3 is similar to the 
elliptic flow and also the mass splitting is observed for 
identified particles. The flat rapidity dependence that re- 
sults from the color flux tubes in the initial conditions is 
in agreement with the observation of the ridge in Ar?-A^> 
correlations. By measuring V3 also in a differential way 
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FIG. 7: (Color online) Pseudorapidity dependence of V3 and 
V2 for charged particles in Au+Au collisions at snn = 200^1 
GeV. 



one might be able to learn something about the amount 
and the size of the initial state fluctuations. 
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